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In a modified ideal Bose gas model we derive an expression for the temperature dependence of 
the superfluid fraction in liquid '*He. This expression leads to a fit formula for the asymptotic 
temperature dependence that reproduces the data significantly better than comparable formulas. 
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I. INTRODUCTION 

The behavior of ■^Hc and ^He liquid on_jDne hand, and of the ideal Fermi and Bose gas on the other hand, strongly 
suggests that there is an intimate relationtl between the A transition in ''He and the Bose-Einstein condensation (BEG) 
of the ideal Bose gas (IBG). Because of the neglect of the interactions one will not expect that the IBG reproduces 
all properties of liquid ''He, in particular not those properties that are directly related to the interactions (like 
specific heat or the compressibility). There are, however, some basic properties of liquid ^He that may be explainec 
by the IBG (like the irrotational superfluid flow). We start by discussing a discrepancy between the IBG and liquid 
^Hc that is — as we will point out — disturbing in view of the suggested intimate relation between the BEG and the 
A transition. 

The critical behavior of the condensate fraction of the IBG is 

^-i^r, /3=J, (1) 

P 2 

where t = [T — T\)/T\ is the relative temperature; the IBG transition temperature is equated with that of th 
transition. The condensate fraction is commonly identifled with the i^uperfluid fraction; this identiflcation explaini 
a number of experimental findings of which the most important oneu is that a superfluid current has no vortices. In 
contrast to Eq. (^) , the experimental superfluid fraction behaves like 

(2) 

P 3 

The suggested intimate connection between the BEG and the A transition is in conflict with fi ^ v. The values 
/? = 1/2 and v k, imply po ^ Ps .just below the transition. 

The standard solution of the conflict between Eqs. (|l|) and (||) appears to be the renormalization-group method. In 
this approach one starts from a Ginzburg-Landau ansatz for the free energy (or enthalpy) that leads to the critical 
exponent 1/2 for the order parameter. For the considered universality class the renormalization procedure yields 
then values near to 1/3 for the critical exponent of the order parameter. This may well serve as an explanation 
of the experimental value of « 1/3 in Eq. (^ but it does not resolve the conflict between Eq. (|^) and Eq. (|^): A 
renormalization is appropriate for the Landau value /? = 1/2 but not for the IBG value /3 = 1/2. The reason is that 
the IBG value is obtained by an exact evaluation of the partition sum. The exact evaluation of the partition sum 
implies a summation over arbitrarily small momenta (or, correspondingly, arbitrarily large distances). Therefore, the 
reasoning behind the renormalization procedure (analytic Ginzburg-Landau ansatz for finite regions or blocks, and 
subsequent transformation to larger and larger blocks) cannot be applied to the IBG free energy. 

Moreover, the critical exponent /3 = 1/2 of the IBG cannot be changed within the IBG frame without destroying the 
mechanism leading to the BEG. The exponent /3 = 1/2 is characteristic for the IBG with the BEG phase transition. 

We have now discussed two points: (i) We expect an intimate relation between the Bose-Einstein condensation and 
the A transition, (ii) The IBG value /3 = 1/2 should be taken seriously (because it is a result of an exact evaluation 
of a partition sum, and because /? 7^ 1/2 is not compatible with the BEG mechanism). From these two points we 
conclude the following: The theoretical /3 = 1/2 in Eq. (^ and the experimental « 1/3 in Eq. (|^) are in conflict with 
each other. 
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Within the frame of the ideal Bose gas model we propose to resolve this conflict by the assumption that noncondensed 
particles move coherently with the condensate. This means that we no longer identify the condensate with the 
superfluid fraction; the condensate is only part of the superfluid phase. A coherent motion can be described by 
multiplying the real single particle functions of noncondensed particles by the complex phase factor of the condensate. 
The superfluid density ps is then made up by the condensate density po plus the density pcoh of the coherently comoving, 
low momentum noncondensed particles. This concept leads to an expression and eventually to a fit formula for the 
temperature dependence of the superfluid density. 

We will stick to the essential characteristic of the IBG (in particular the EEC) but introduce some modifications 
(for example, Jastrow factors) that are necessary for a realistic approach to liquid *He. This modified IBG is called 
almost ideal Bose gas model (AIBG). The AIBG has been introduced some years ago as an attempt to explain the 
(nearly) logarithmic singularity of the specific heatu. Some consequences of the decomposition Ps = Po + Pcoh have 
been discussed in Rcfs. |^ and ^. The present paper is devoted to the investigation of the temperature dependence of 
the superfluid density in this model. The necessary details of the underlying model, the AIBG, will be given below. 

The form of the temperature dependence of the superfluid fraction is derived in Sec.|l[ This leads to a fit formula for 
the temperature dependence of ps that is applied to experimental data and compared to other fit formulas (Sec. III). 
Sec.0 discusses the temperature dependence of the condensate density. Section ^ presents scaling arguments on 
the basis of an effective Ginzburg-Landau model; this includes a qualitative explanation of the coherent comotion of 
noncondensed particles and leads to restrictions for some of the parameters of the fit formula. 



II. AIBG FORM OF THE SUPERFLUID FRACTION 



A. Many-body wave function 

Following Chesteifl we multiply the IBG wave function \I'ibg by Jastrow factors = 11 /u > 

N 

i<j 

We consider i = 1,2, N atoms in a volume V. The occupation numbers are parameters of the wave function; in 
physical quantities they are eventually replaced by their statistical expectation values (nk)- The Jastrow factors take 
into account the most important effects of the realistic interactions; with, a suitable choice for the fij (for example, 
fij{r) = exp[— (a/r)^] with a and b determined by a variational principle^) the wave function (H) leads to a realistic 
pair-correlation function. 

The IBG wave function ^'ibg in Eq. (^) is the symmetrized product of single-particle functions. We display this 
structure admitting at the same time a phase field $ of the condensate: 

vl/ = 5i^[exp(i$)]"'' n [v^kT . (4) 

Here S denotes the symmetrization operator. The c/Jk are the real single-particle functions of the noncondensed 
particles. The schematic notation [t/Sk]"'' stands for the product (/7k(ri/+i)-</?k(ri/-i-2)'- ■ ■■^ki^v+nu.)'-, this notation applies 
also to [exp( i $)]"<' . All no condensed particles adopt the same phase factor exp( i <I'(r)) forming the macroscopic wave 
function 

^(r) = y^exp[ici>(r)] . (5) 

The phase field, 4> describes the coherent motion of the condensate particles. (Actually, one has to construct a suitable 
coherent statecl. This point is, however, not essential for the following discussion.) This motion is superfluid if the 
velocity Ug = ftV^/m is sufficiently small. r, 

Equations (^ and (||) are a well-known descriptiona for a superfluid motion in the IBG. In this description the 
superfluid fraction ps/p equals the condensate fraction hq/N = po/p. The role of the Jastrow factors in this context 
will be discussed in Sec. |IIB| . 

In order to dissolve the discrepancy between Eqs. (0) and (^) we assume that noncondensed particles move coherently 
with the condensate. This is possible if noncondensed particles adopt the macroscopic phase of the condensate: 

*=5^^[exp(ici>)]"'' Y[ [^kexp(i $)]"'= J] [^k]"" • (6) 

< fe < fccoh k > fccoh 
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We assume the phase ordering for all states with momenta below a certain coherence limit fccoh- For the low lying 
states with ^ 1 such phase ordering is relatively easy beca use i t requires only a small entropy decrease. At this 
stage, kcoh should be considered as a model parameter. In Sec. V A the existence and the size of this coherence limit 
will be made plausible. 

We evaluate particle current for the wave function (||): 



j,(r,nfe) = ( 



N 



+ c.c. 



p h 
N m 



no 



(7) 



The prime at the sum over the momenta k < fccoh means that the k = contribution is excluded. In coordinate space, 
the current operator reads j„ = — i/iV„/(2m) + c.c. It acts on all ri-dependences. Because of the added conjugate 
complex term all contributions from the real functions (the in F or the (pk) cancel. The only surviving terms are 
those where j„ acts on the phase $. 

For a superfluid motion with Ug — hV^/m and in the statistical average, js of Eq. equals ps Ug. We may then 
read off the superfluid fraction. 



P 



N 



PO + Pcoh 
P 



(8) 



This expression will be evaluated in Sec. II C 

The ansatz (|^) leading to Eq. shows in which way noncondensed particles may contribute to the superfluid 
density. 



B. Condensate density 



We discuss in some detail what is meant by the terminus "condensate density" , in particular with respect to the 
Jastrow factors in Eqs. (^) or (||). 

The exact condensate density may be defined by 

(vl/|0+(r)0(r')|*> po-ct ^ (9) 

where is the exact many-body state and the (j)'^ and 4> are single-particle creation and annihilation operators. For 
finite temperatures one has to take the statistical expectation value of pQ^'^^ (we do not introduce a different symbol). 

For an IBG wave function V&ibg the condensate density is given by p™°'^''' = no/V. In the statistical average this 
this model condensate density becomes 

model ("-o) /-.^.x 

Po = — ■ (10) 

The exact many-body state in Eq. (^ may be approximated by Eq. (H), or by ^' w -F for the ground state. -Ln-this 
case the relation between both condensate densities is well known: The model condensate fraction is depletedt3't3 by 
the Jastrow factors F, for example, from p[f°'^°Vp = 1 to Po™'^*//c « 0.1 for T = 0. 

The above calculation leading to Eq. (||) demonstrates the following point: In contrast to the density the 
current density p^°'^^^ Ug is not depleted. The reason is that in Eq. (|^) all derivatives of the real Jastrow factors cancel 
(because of the added conjugate complex term). 

On the basis of this point we arrive at the following statements about the role of the densities Pq^^^*^, and 

Ps- 

(a) Since the current density pg""'*'^' Ug is not depleted we may identify (and not Po''^'^*) with the square j^/^p 
of the macroscopic wave function. Irrespective of the Jastrow factors we may use Eq. (^ as it stands. For 
a superfluid flow, the phase $(r) of the macroscopic wave function (^) fixes the velocity field Ug = hV^/m. 
The basic relations for the superfluidity (like curlug = and the Feynman-Onsager quantization rule) are not 
affected by the Jastrow factors in the many-body wave function. 

(b) The exact condensate density is a quantity of its own right. It is the density of the zero momentum particles in 
the liquid helium. Recently Wyattt3 reported about a rather clear experimental evidence for this condensate. 
For a review about the attempts to determine pg'^^'^* experimentally we refer to SokolEj. 
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(c) The assumption that noncondensed particles move coherently with the condensate is introduced by the step from 
Eq. (^) to Eq. (^). Again, this step does not alter the basic relations following from Eq. (^) (like Ug = ft V<l>/m, 
curlug = and the Feynman-Onsager quantization rule). 

(d) For T = the value Po"°'^°Vp — 1 yields ps/ p = 1 [as we will see, pcoh in Eq. (H) contributes only in the vicinity 
of T\\. In contrast to this the connection between p'^^^^ / p ~ 0.1 with ps/ p = 1 is less obvious. For T w the 
value P™°'^°Vp ~ 1 implies Ps/ p ~ 1. For describing 1 — Ps/ p quantitatively one must however include phonons. 
This is not done in Eqs. (Q) or ^ because our primary object is the asymptotic temperature region. 

We summarize this subsection: As far as the superfluid current is concerned the model condensate density is not 
depleted. The model condensate density is the fundamental constituent of the superfluid density. 

In the following the model condensate density will again be denoted by po and called condensate density. 



C. Superfluid density 

We evaluate the expression (^) for the superfluid density. Our model assumes expectation values {rik) that are of 
the IBG form, 

^"'^ " exp[(efc-Ai)/fcBT]-l " exp(x2+T2)-l " ^^^^ 

Here p is the chemical potential, = h?k'^/2m are the single-particle energies, and /cb is Boltzmann's constant. In 
the last expression we introduced the dimensionless quantities = —^./k^T and 

A|k| . , , 27rft 
X = — , with A = -^=^= . (12) 
\/47r v27rmfcBT 

The transition temperature of the IBG is given by the following condition for the thermal wave length A: 

A(T,)=[<(3/2)]'/^ (13) 

where C(3/2) = 2.6124 denotes Riemann's zeta function. In applying our almost ideal Bose gas model (AIBG) to the 
real system we identify T\ with the actual transition temperature. In the following we use the relative temperature 

t=^. (14) 



We evaluate the condensate density: 



p ^ N ' C(3/2) ^ ' 

Riemann's generalized zeta function is given by gp(r) = exp(— nr^)/?!-^, and C(p) — gp{0)- The chemical potential 
p or, equivalently, r may be expanded for \t\ <^ 1: 



_ ( at + bt^ + ... it>0) 
^^*>~\lk^-\a'\t\+b't^ + ... it<0) ■ ^^^> 

For t > Eq. ([l5|) with po/p — yields (1 + t)^^'^ g3/2{''') — C(3/2). This condition determines the temperature 
dependence of T{t) and in particular the coefficients a, b, . . . , for example a — 3 C(3/2)/(47r^/2). 

For t < the IBG yields t = 0. In the AIBG we admit nonvanishing coefficients a', 6',. . . in Eq. (p^. This makes 
the expansion (|l^) more symmetric; it corresponds to a phenonemological gap between the condensate level and the 
noncondensed particles. A coefficient a' ^ does notjaffect the BEG as the most important feature of IBG. It avoids, 
however, the divergence of the static structure factortj S{k) for fc ^ and greatly improves the unrealistic (oc T'^/^) 
behavior of the specific heat. In view of the successful roton picture it is not too surprising that a gap is necessary 
for a quantitative description of the superfluid density (or of the specific heat). As we will see, a realistic description 
of liquid helium requires a' ~ 3; the next coefficient b' will not be needed. 

From Eq. (M) and with Eq. (Ill) we obtain 
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p 

where 



^° + (i<0) (17) 



We evaluate now the density of the comoving particles 



Pcoh _ ^/ (rifc) _ 4(1+^)3/2 x-^dx 



p iV 0FC(3/2) 7o exp(a;2 + r2)-l 



Using y/[exp(y) — 1] = 1 — y/2 + y /12 ^ ... we obtain 

ftoh _'l(l + i)"V. , a:coh i=„h 



The convergence of this expression is excellent; for the actual parameter values and for \t\ < 0.1 the terms not shown 
are of the order 10~^. 

We have not yet specified the coherence limit fccoh- For \t\ <C 1 we will find ps ^ Pcoh ^ fccoh for the superfluid 
density and ps fc^oh ~ ^coh ^^r the kinetic energy of the fluctuations. Requiring that this kinetic energy scales with 
the free energy F ~ — (no)^ ~ — yields 

fccoh ~ . (21) 

This scaling argument will be presented in more detail in Sec. 0. 

Inserting Eq. (^l|) in Eq. ([20|) and using Eq. (17), the superfluid fraction contains the powers \t\, and so 



P P 



D. AIBG assumptions 

We summarize in which points the AIBG, the almost ideal Bose gas model, deviates from the IBG: 
1. The IBG wave function ^'ibg is multiplied by Jastrow factors, 5" = F^'ibg- This is a well-known approac 



J. 



2. By the symmetric expansion ( p^ we admit a gap between the condensed and the noncondensed particles. This 
modification preserves the most basic features of the IBG, in particular the BEG mechanism and the critical 
exponent /? = 1/2. 

3. The noncondensed single-particle states below the coherence limit fccoh adopt the macroscopic phase of the 
condensate. The leading exponent for the coherence limit fccoh is determined from a scaling argument. 

III. FIT TO EXPERIMENTAL DATA 

We compare the temperature dependence of our model expression for the superfluid density with experimental data. 
The model expression contains unknown parameters; it provides a fit formula for the data. It will turn out that this 
fit formula is significantly better than comparable fit formulas. 
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FIG. 1. The asymptotic model fit (MF), (Ps/p)mf ~ a,i + 0,2 |i| + "3 1^1*''^, reproduces perfectly the experimental 

superfluid fraction (ps/p)cxp- The solid points mark the relative deviations between the fit formula and the experimental values 
of Ref. |l^. As a guiding line we draw straight lines between neighboring points. The broken line corresponds to two standard 
deviations. 



A. Asymptotic temperature range 

1973 data by Greywall and Ahlers 



A restriction to the first three terms in the expansion (22) yields a three-parameter model fit (MF) 

Pl^a^\t\V^ + a2\t\+a:,\t\''^ (MF) . 
P 



(23) 



Figure |i| shows that the MF yields an excellent reproduction of the data by Greywall and Ahlersll3 for saturated vapor 
pressure. We used all data points with temperatures \t\ < 0.03. For a minor improvement of the fit we shifted the 
temperature values by —0.5 x 10~^; this is well below the experimental uncertainty oi 5t = 2 x 10~^. The parameters 
of the fit shown in Fig.|l] are: 



ai= 2.3233, 02 = 1.0258, 
As an alternative we consider the standard fit (SF) 



03 = -2.0065 . 



(24) 



D\t\ 



(SF) 



(25) 



P 

which is used by GrjeMisiall and AhlcrsEl, and that is motivated by the renormalization-group theory. The fourth 
parameter A is oftenll3ll3 set equal to 1/2 because the fit is not very sensitive to it. We will use the SF with A = 0.5 
as a three-parameter ansatz. 

The fit parameters are found by minimizing the sum of the quadratic deviations, 



Ps 

PJht 



(Ps/P)fit 
(Ps/P) 



cxp 



(26) 



Here A^d is the number of data points and a^ei is the relative standard deviation. The standard deviation a for 
(Ps/p)fit ~ (Ps/p)oxp is given hy W ~ 1/ct^. The dominant experimental error is that in the temperature. This leads to 
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the weightll3 W = |fp/^/(5|tl^, where 5\t\ = max (2 x 10"^, 10~^ |t|) is the temperature uncertainty. The corresponding 
2 tJi-ei hne is shown in Fig.jlJ. 

In the given form both MF and SF (with A = 0.5) are three-parameter fits. We compare both fits by calculating 
their ratio: 

w 8.8 (data for \t\ < 0.03) . (27) 

Xmf 

As seen from Fig.|^ the MF reproduces the experimental data {x^/N^ ~ 1.10). The large ratio ( |27| ) means that the 
SF does not reproduce the data in the considered temperature range. 

We remark that the SF fits the data in the considerably smaller range |t| < 0.004. This smaller range is used in Ref. 
[isl , presumably because it was realized that the SF does not fit the data in the larger range. For a three-parameter 
fit the range |t| < 0.004 appears to be rather small; we note that already a one-parameter fit (aijip/'^) reproduces the 
data within 2% in the relatively large range |t| < 0.08. |— ■ 

We considered also the data at higher pressures by Greywall and Ahleralj. Here we found ratios xIf/Xmf between 
1 and 2, and values of XMpZ-^d in the range between 3.6 and 15. This means that the MF is only slightly better than 
the SF without yielding satisfactory fits. This is (at least partly) caused by jumps in the experimental data points. 
For example, compared to a smooth fit curve (SF or MF or any reasonable fit formula) there is a jump of more than 
ten standard deviations between the data points (|t|, pjp) (0.0014391, 0.028 144) and (0.0012631, 0.025 624) for 
P = 7.27 bar. 



1993 data by Goldner, Mulders and Ahlers 

Newer measurements of the superfluid density are reported by Goldner et al.lli and by Marek et al0. We consider 
the data byj-Qoldner et al. because these authors published an explicit data list. 

The datallj extend to about \t\ = 0.01; all these data are used for the fits. Fig.H shows how the three-parameter 
MF reproduces these data. We discuss this result in a number of points: 

1. Obviously the scatter of the data is generally larger than the estimated error (listed as 5ps/p in Ref. |l^, and 
called fj in Fig.^. There are several jumps of the size of ten standard deviations; the most dominant jump 
(between the values for \t\ — 0.00031910 and \t\ = 0.00039793) is about 30 times larger than the estimated error. 
This statement is basically independent of the fit formula used (see also Fig.^). It is extremely unlikely that 
the actual superfluid fraction contains such jumps. The different sizes of the jumps restrict the possibility to 
discriminate between various fit formulas. This is also the reason why we considered first the older 1973 data 
by Greywall and Ahlers. 



2. The three-parameter SF yields a significantly larger x^ value: 



XSF 



2.7 



(28) 



3. Goldner et alJll used the following extended standard fit (ESF) 

^ = ko\tmi + D\t\^)il + h\t\) (ESF) 



(29) 



4. 



with A = 1/2. Using the same parameters as in Ref. we obtained Y^gp/xMP ~ l-^- This might appear as a 
small difference between MF and ESF. A comparison between Figs. ^ (MF) andJ3| (ESF) shows, however, that 
the MF does a better job although it has one parameter less. Goldner et al.ll3noted that there is a serious 
discrepancy between the ESF and the data, in particular in the range \t\ ~ 10~^ to 10^^ (their Fig. 17). The 
comparison between Fig.|| and ^ shows that this discrepancy is significantly smaller for our model fit. This 
improvement is not so evident in the x^ ratio because the x^ values are on a high level for any fit formula (due 
to the jumps). 

Looking at the scatter of the data one might tentatively assume a standard deviation that is five times larger 
than the one assumed. Drawing then a new 2 a line the discrepancies in the range \t\ = 10~^ to 10~^ in Fig.^ 
may be judged as not very significant. They may, however, hint at an unexplained structure in the temperature 
dependence of the superfluid fraction. 
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FIG. 2. The asymptotic model fit (MF), Eq. (||), is applied to the 1993 data of Ref. The deviations between the fit 
and the data are given in units of the experimental error; the broken line corresponds to two standard deviations. For some 
temperatures there are several experimental ps values (each leading to a full point); in this case the guiding line goes through the 
average value of the deviation. Sometimes the experimental temperature values are close together (but different); this results 
in apparently vertical pieces of the guiding line. From the figure it is obvious that (i) the statistical errors are considerably 
larger than the assumed a, and (ii) the data contain several large jumps between neighboring points. 
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FIG. 3. The extended standard fit (ESF), Eq. (|9|), is applied to the 1993 data of Ref. |l|. The presentation is the same as 
in Fig.H The systematic deviations between the fit and the data are significantly larger for this four-parameter ESF than for 
the three-parameter MF (Fig.^). 
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B. Extension to lower temperatures 



We apply the model expression for the superfluid fraction in the temperature range 1.2K < T < Ta, where \t\ is no 
longer much smaller than 1. For this purpose we use the model expressions ( |l5|) and (jl^) for po and pcoh, respectively, 
and expand t and Xcoh (rather than ps itself) into the relevant powers of \t\. 

The expansion t — a' \ t\ + b' t'^ + . . . may be broken off after the first term because t ^ corresponds to a gap and 
leads to an exponential decrease [oc exp(— r^)] of the noncondensed contribution in Eqs. (|l5|) and (LO]). Therefore, the 
noncondensed contributions become rather small before the next terms in the expansion for r contribute significantly. 

As far as the coherence limit fccoh is concerned we have no information about the continuation of Eq. (^) into an 
expansion. In view of the success of Eq. ( |2^ ) we will certainly not admit exponents that would violate the form ([2^). 
In accordance with Eq. ( p^ we may admit the form Xcoh = xi + X2 |i| + ^3 \t\'^^^ + . . .. This expansion may be 
broken off, too, because Pcoh of Eq. ( p!9| ) is damped exponentially [cx exp(— r-^)] for increasing \t\. Including the terms 
with the parameters xi, X2, and X3 preserves the variability for the parameters ai, 02, and 03 in Eq. (|2^). 

Due to the exponential damping of the noncondensed contributions a cut in the expansions for r and a:coh leads 
much further than a cut in the expansion for ps/p itself. In this way we arrive at the following unified model fit (UMF) 
formula: 

p C(3/2) V^C(3/2) Jo exp(x^+r^)-l 

with 

T^a'\t\, Xcoh = niax(0, xi|i|2/3+x2|t|+X3|i|4/3) . (31) 

As we will see, this formula provides a unified description of the asymptotic region as well as of the less asymptotic 
(the "roton") region. 

The parameters xi, X2, and X3 are related to the oi, 02, and 03 in Eq. (|2^) and essentially fixed by the asymptotic 
region. We have restricted a;coh explicitly to non-negative values because the expression xi \t\'^^'^+X2 \t\+X3 might 
become negative for larger \t\ values [where, however, the density pcoh tends to zero anyway because the exponential 
decrease cx exp(— r^); see also Fig.||]. l— . 

For a fit in the ra»ge 1.2K < T < Tx we combined the data by Greywall and Ahlerslia for \t\ < 0.04 and that 
by Clow and ReppyEJ (run IV) for \t\ > 0.04. At \t\ — 0.04 both data sets are compatible with each other. The 
systematic errors and the deviation due to slightly different prsssures (roughly 1% between saturated vapor or normal 
pressure) just happen to cancel each other. Clow and Repp}£3 remark that their "values of Ps/p have a scatter of 
about 1/2%"; we interpreted this as (Trci — 0.005 for our fit [i.e. for the minimization of Eq. (Pq)]. 




0.0 0.1 0.2 0.3 0.4 \t\ 

FIG. 4. The unified model fit (UMF), Eq. (|^), is used to reproduce the experimental superfluid fraction (ps/p)cxp in the 
temperature range from Tx to 1.2 K. For \t\ < 0.04 we use the Greywall- Ahlers data (Ref. |l^, for jtj > 0.04 that of Clow and 
Reppy (Ref. ^8|). The figure depicts the deviations for the range \t\ > 0.04; for \t\ < 0.04 the deviations are quite similar to 
that shown in Fig. hi The broken line corresponds to two standard deviations. 
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A fit of the combined datallllli leads to a result that is quite similar to Fig. |] for |t| < 0.04 and that is shown in 
Fig. ^ for \t\ > 0.04. The fit parameters are: 

a' ^ 3.0380 , xi = 2.6998 , X2 = -0.8063 , = -3.9631 (UMF). (32) 

Alternatively we may use the parameter /, Eq. (p^, and calculate the parameters ai, a2, and following from the 
asymptotic expansion of Eq. ( |30| ) : 

/ = 5.6225, ai = 2.3323, as = 0.8035 , 03 = -0.4704 . (33) 



If an expansion is broken off as in Eq. ( p3[ ) the last term tries effectively to simulate the missing terms. Since the 
UMF supplies higher-order terms it is not surprising tiiat the last coefficients in Eqs. (|4|) and (|3^) are quite different. 

Alternatively we used the data by Tam and AhlersO that extend, however, only down to 1.5 K. This yields similar 
parameter values. 

The standard fit for temperatures above 1 K but excluding the asymptotic region is the two-parameter roton fit 
(RF), 

^ oxpf-A^) (RF). (34) 



P Vt \ kBT 
Using the data of Ref. |l^ (run IV) we obtain 



XuMF 



for 1.2K < T < 2.07 K. (35) 



This ratio is reduced to 2 if we restrict the temperature by T < 2K. These ratios imply that the unified model 
expression is quite good for intermediate temperatures, too. 

The RF is based on Landau's quasiparticle model that cannot be extended to T\ without loosing its physical basis. 
The standard description for the range 1.2 K < T < T\ would be a combination of the SF ( p5| ) and the RF (jsj). In 
contrast to this, our model provides a unified fit (^) in this range. Although containing one parameter less (than the 
combination of SF and RF) this unified fit is superior to the standard description. 

As already mentioned, the expansion ( |l6|) implies a gap between the condensed and noncondensed particles. This 
gap appears to be essential for the reproduction of the data in the intermediate range T > 1 K. This gap should in 
some way be related to the roton gap A. This relation cannot be expected to be simple and obvious because one gap 
belongs to a model (Landau) for T <^T\ and the other to a model (AIBG) for T ^ T\. We note that our gap vanishes 
for T ^ T\, and that the roton concept becomes less sharp for increasing temperature (for T = IK the widths of 
roton states are already comparable to their energies). 

For T <^T\ Landau's quasiparticle model is, of course, the right model. The model fit (|o|) yields still reasonable 
values for Ps/ p but it must fail in the quantitative reproduction of 1 — Ps/ p because the phonons are not described by 
the wave function (||). 

IV. CONDENSATE DENSITY 

The unified model fit, Eq. ( |30|) with Eq. (^^, defines the decomposition of the superfluid density into the condensate 
density and the coherently comoving density. The temperature dependence of this decomposition is displayed in Fig.||. 
In this section we discuss in particular the temperature dependence of the condensate density. 

The contribution of pcoh is decisive near T\ but negligible for lower temperatures. The comoving density pcoh carries 
some entropy because it does not correspond to a single quantum state. This entropy content is quite small because 
it is due to the lowest single-particle states with (rxk) ^ 1. It is below the present experimental limits but should be 
detectable; for thes e poin ts we refer to Refs. || and ||. 

As shown in Sec. Ill B , the expression t = a' \t\ works quite well for fitting the data down to about 1.2K. Inserting 
T = a' |i| in Eq. (|l5|) yields 

ei^l-^, + tfn?2lg^. (36) 
P C(3/2) 

Using a' of Eq. (|3^), this temperature dependence is shown by the dashed line in Fig.||. 
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FIG. 5. Decomposition of the superfluid density ps into the (model) condensate density po and the coherently comoving 
density pcoh as a function of the temperature T. For T T\ the comoving density pcoh oc \t\'^^^ is the dominant contribution 
to the superfluid density. The full line represents the model expression for ps. In the given scale, this curve coincides with the 
experimental superfluid density for all temperature values. 



The asymptotic expansion of Eq. (B(|) reads p^/ p ^ f\t\, Eq. (17), where 



The numerical value is taken from Eq. 



3 

We found that 



P 




(37) 



(38) 



may be used as an approximation for Eq. (|3^). The maximum relative difference between Eqs. ( |38D and ( |36D is about 
2%. For T ^Tx both expressions, and (|6|), yield po/p ~ f\t\. 

The right-hand side of Eq. (^) is an old fit formula for the superfluid fraction Ps/p (for example, Fig. 27 of Ref. |^). 
In the framework of our model, this historic fit formula may be interpreted as the approximation ps ~ po ■ The obvious 
shortcomings of Eq. ( p^ ) as an approximation for Ps/ p are the following: (i) The neglect pcoh leads to a qualitatively 
wrong asymptotic behavior (difference between the full and the dashed line in Fig. ^) . (ii) The step from Eq. ( |3^ ) 
to Eq. ( p8| ) as well as the use of r = a'\t\ make the expression an approximate one already for po/p. (iii) For small 
temperatures 1 — Ps/ p = {T/TxY is quantitatively wrong (because the phonons have not been taken into account). 

We consider once more the exact condensate fraction p™^'^*/p introduced in Sec. [IB. We denote its value at T = 
by Uc- Assuming that the depletion of the condensate (from 1 to ric ~ 0.1) is temperature independent we obtain 



„cxact 



PO 
P 



1 



Tf 



f 



(39) 



as an approximate expression for the temperature dep| 
temperature dependence is given in Fig. 2 of Snow et al.t 
expression (^) agrees with the data. 



idence of the exact condensate fraction. The experimental 
J. Within the relatively large experimental uncertainties the 



V. EFFECTIVE GINZBURG-LANDAU MODEL 



In our approach, the macroscopic wave function 



V'(r) = Jy [i = \/po exp [i$(r)] 



(40) 



plays the role of the order parameter. We investigate the free energy as a function of this order parameter. 
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A. Coherence limit 



We start by presenting a qualitative argument for the existence and the meaning of the coherence Umit kcoh- 
The macroscopic wave function ^ may contain equiUbrium and nonequihbrium excitations. A superfluid motion 
with Us = (/i/m)V$ is a nonequihbrium excitation. At finite temperatures, there are thermal fluctuations of the 
order parameter, i.e., equilibrium excitations. We consider the average momentum of these fluctuations. 



fcfluct = |V$| . (41) 

The bar denotes the statistical average. The momentum /cfluct will be a function of the temperature. It is related to 
the correlation length ^ « 1/fcfluct- 

The single-particle states are described by real functions tpk in Eq. (||). We consider the possibility of phase 
fluctuations for a low-lying state with ^ 1, too. After the replacement (/Sk — > <i5kexp(i$k) in Eq. (||) these 
fluctuations may be described by the fields <i>k(r). 

Let us first assume that the additional phases vanish, <i>k ~ 0. In this case, the average kinetic energy ft^fc^uct/^"^ 
of a condensed particle would exceed that of a noncondensed particle with k < fcfiuct- The energy sequence of the 
single-particle states is, however, a prerequisite of the BEC; the condensate must be formed by the particles with the 
lowest energy. In order to preserve the energy sequence of the low-lying states we require the phase ordering 

$k(r) = $(r) for fc < • (42) 

This argument does not apply to the states with higher momenta. By this qualitative argument we obtain the 
many-body wave function (^) with 

^coh = fcfluct • (43) 



B. Free energy 



The statistical expectation value pq ~ \t\ can be obtained by minimizing the common Landau energy Fi^/V = 
Rt\xi}'\^ -I- (with regular coefficients R and U). The fluctuation term Ffluct/^^ = (fi^/2m) |V-0P equals the 

kinetic energy density p^v? jl of the condensate only; here u = {h/m^jV^. The phase coherence assumed in Eq. (^) 
implies that pau^ must be replaced by PsU^. This leads to the following effective Ginzburg-Landau ansatz 

^GL ^ WKFl ^^P_^ , +Rt\M'' + U ■ (44) 
V V 2m po 

Assuming that the leading exponent of Xcoh is not greater than 1, Eq. (|^ yields 

Pcoh ~ a;coh ■ (45) 

The equilibrium fluctuation term becomes then 

i^fluct OC Ps fcfluet ^ iPO + Pcoh) fccoh ^ PO X^oh + ^coh ■ (46) 

The asymptotic form of the Landau part of the free energy behaves like 

Fi^(xRtpo + Up^ (47) 

We require now scaling invariance. This means that i^fluct must have the same leading \t\ dependence as Fi^, i.e., 
Po^coh + ^coh ~ From Pox^^^^ ~ we would obtain Xcoh \t\^^^ and x^^^^ ^ \t\'^^^ in contradiction to the scahng 
assumption. Therefore, scaling requires x^^^i ~ I^P o^' 

fccoh - . (48) 

This implies ps Pcoh ^ for the superfluid density and ^ l/fcfiuct ~ fni:.the correlation length. 

The mass coefficient Ps/po ~ 1^1"^^'^ in Eq. (^J) is singular. Ginzburg and SobyaninEj have introduced a comparable 
effective Ginzburg-Landau model with nonanalytic coefhcients, too. In Ref. ^the nonanalytic coefhcients (like \t\'^^'^ 
for the IV'P term) are phenomenologically introduced in order to reproduce the right critical exponents. 

The divergent mass coefficient Ps/po ^ \t\~^^^ damps the critical fluctuations such that Eq. ( ^ ) becomes scaling 
invariant. In this sense, the model ( ^ ) has properties similar to the common Ginzburg-Landau ansatz in d = 4 
dimensions. This means that Eq. (|4^) might be used down to \t\ = and that the critical exponent of ps might be 
indeed exactly 2/3. This possibility is supported by the excellent fit obtained for Eq. (p^). 
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C. Further scaling restrictions 



The equilibrium Landau free energy contains integer powers of t only: 

Fi^ ^ ...t^ + ...t^ + ...t"^ + ... . (49) 

The asymptotic form of the superfluid density ( p^ is compatible with the expansion Xcoh oa Xi \t\'^^^ + X2 \t\+X3 + 
. . .. In the fluctuation term this expansion will, however, in general lead to noninteger exponents: 

i^fluct CX Ps fccoh = (PO + Pcoh) fccoh - -t^ + - l^r^' + - 1^1'^' + - + - • (50) 

Scaling for Eq. ( ^ ) implies also that the amplitudes of the nonanalytic terms vanish. This condition yields relations 
between the expansion parameters a', 6', . ■ ■ Sind xi, X2, X3, . . . that may also be expressed by the coefflcients in 
Eq. ([2^). The condition of a vanishing amplitude of the li]^/"^ term can be evaluated straightforwardly and yields 

^-0.58 or a2-l. (51) 

8 

These theoretical values compare well with the fitted values given in Eqs. ( p^ or (|3^). 
The condition of a vanishing amplitude of the \t\^/^ term yields 

.3 = M!-|:i ,52, 

04x1 3xi 

and a corresponding expression for a^. These relations are not fulfilled by the parameter values found in the fits. 
The reason is probably the following: The expansions (|2^) and (31) are cut after the term. In a fit it is then in 



particular the last term that tries to simulate the neglected terms. 

VI. CONCLUDING REMARKS 

We have modified the IBG in such a way that it might be applied to liquid helium. We summarize the novel views 
and main results of our approach. 

1. The IBG value (3=1/2 for the critical exponent of the condensate should be taken seriously. It is not subject 
to renormalization because it results from a calculation that already includes a summation over arbitrarily large 
lengths, and it is essential for the BEG mechanism. 

2. The model condensate contributes fully to the superfluid density; it is not depleted by the Jastrow factors. 

3. In order to reproduce the critical exponent « 1/3 of the superfluid density we have assumed that noncondensed 
particles below a certain momentum k^oh move coherently with the condensate. The coherence limit fccoh has 



been made plausible in Sec. VA 



The contribution of noncondensed particles to the superfluid density offers a solution of the so-called macroscopic 
problen£3 of liquid helium. This problem reads as follows: If the superfluid density corresponds to single 
quantum state {ps oc IV'P) then the approach to an equilibrium state [with ps = Ps{T)] cannot be understood. 

4. We have derived a fit formula for the temperature dependence of the superfluid density. This fit formula 
reproduces the data significantly better than comparable expressions. This feature as well as qualitative scaling 
arguments suggest that the critical exponent v of the superfiuid density might be exactly equal to 2/3. 

5. The temperature dependence of the decomposition of superfiuid density into the model condensate density and 
the coherently comoving density is given. A simple formula for the temperature dependence of the depleted 
condensate density is presented. 
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